Effect of external conditions on the structure of scrolled graphene edges 
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Characteristic dimensions of carbon nanoscrolls — "buckyrolls" — are calculated by analyzing the 
competition between elastic, van der Waals, and electrostatic energies for representative models of 
suspended and substrate-deposited graphene samples. The results are consistent with both atomistic 
simulations and experimental observations of scrolled graphene edges. Electrostatic control of the 
wrapping is shown to be practically feasible and its possible device applications are indicated. 



Carbon nanoscrolls (CNSs) are intriguing materials 
that differ from both carbon nanotubes and carbon 
nanoribbons. Unlike nanotubes, CNSs have no caps nor 
their electronic states are subject to periodic boundary 
conditions. Unlike nanoribbons, CNSs are curved, and so 
they respond differently to uniform external fields. CNSs 
have additional mechanical degrees of freedom, e.g., the 
inner and outer diameters whose manipulation may be 
utilized in many new applications. 

Fabrication methods of CNSs include arc discharge, 1 
sonication of graphite, 2 ^ scratching it with an AFM tippJ 
electrodeposition of graphene in a gaseous atmosphere, 5 
and immersion of graphene in alcoholP Theoretical stud- 
i e ^Ztt!4] indicate that the lowest energy configuration of a 
large graphene flake is fully wrapped. However, it is sep- 
arated from a fiat state by a large energy barrier due to 
bending rigidity. This is one reason why flat free-standing 
graphene membranes can exist In this paper we inves- 
tigate another possible reason: the wrapping is arrested 
by unfavorable boundary conditions or external fields, 
and so what is observed is a partially scrolled edge of the 
otherwise mostly fiat sample, see Fig. [I] We estimate 
dimensions of such edge CNSs, study the charge distri- 
bution therein, and discuss how they can be controlled 
electrically. Our approach relies on the continuum elas- 
ticity and classical electrostatics theories. It should be 
adequate for typical CNS, which contain anywhere from 
several tens to many thousands atoms around the cir- 
cumference. 

We begin with revisiting the problem of a CNS of a uni- 
form cross-section with no axial elastic stress. The main 
structural property of interest is the smallest possible in- 
ner diameter d m i n of a stable CNS. The input i nto o ur 
calculation are the following material parameters ISllZHUD 
which are derived from the properties of graphite and car- 
bon nanotubes: the interlayer distance h = 0.34 nm, the 
elastic bending modulus D = 1.5 eV, the two-dimensional 
(2D) Young's modulus Y 2 = 340 N/m w 2100eV/nm 2 , 
and the van der Waals energy per unit area of contact 
7 = (40 ± 10) meV/atom = (1.5 ± 0.4) eV/nm 2 P 

The scrolling is driven by the competition between van 
der Waals attraction of graphene to itself and its bending 
rigidity. For a large CNS the former dominates and the 



scrolling is advantageous. We find that the difference in 
energy e per unit length of a CNS and a fiat sheet of the 
same area area has an approximately linear dependence 

e{s) « -(s - s )g, s > s , (1) 

on the total arc length s. In the practically important 
range 10 < s (nm) < 50 we get g = 1.2eV/nm and 
so = 12 nm. The coefficient g is somewhat smaller than 
7 because the bending energy is still significant at such 
s. At s ^> so, the derivative de/ds tends to —7. 

To obtain these results we assume that the CNS is 
shaped as an Archimedian spiraP^ = a0, where 
< <j) < is the polar angle in the cross-sectional plane 
and a = h/(2it). Ignoring any changes in the electron 
spectrum, we take e = e v dw + £b where sb and e v dw 
are the bending and the van der Waals energy, respec- 
tively. Expressing those in terms of the local curvature 
and the appropriate arc lengths, we get the energy func- 
tional e(Q,&) which can be easily minimized numerically. 




FIG. 1. (Color online) Left: a fragment of Fig. SI from Ref.[T6l 
showing an edge CNS, on which we overlaid the x-y axes to 
bring out its parabolic shape. Bottom: an illustrative sketch 
of the same structure. Right: another suspended graphene 
device with scrolled edges (courtesy of A. K. Geim). 
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We will now calculate T and use it to show that 



FIG. 2. (Color online) Energy and the inner diameter d — 
2r + (h/2) of a uniform CNS computed within the continuum 
elasticity model (thick lines). The markers correspond to the 
smallest stable CNS. The straight thin line is Eq. 



The results are shown in Fig. [2j Thus, the arc length 
Smin = 12 nm of the smallest stable CNS is determined 
by £(s m i n ) = 0. The corresponding inner diameter is 
^min = 2.2 nm, see Fig. |2j Our calculation avoids several 
approxim ations m ade in previous elasticity theory mod- 
els of CNS 8 13 21 and as a result gives somewhat different 
d m [ n for the same input parameters. With the quoted 
values we ach ieve a nominal agreement^ with molecular 
dynamic J^U and ab initio simulations P^l 

As seen in Fig.[2j there is a large energy barrier e(0) ~ 
5 eV/nm that separates the initial flat and the final curled 
states. This is why the optimal inner radius can hardly be 
achieved in practice.^ Realistic r ^> d m [ n /2 are presum- 
ably determined by interaction of graphene with scrolling 
"chaperones," e.g., stray fluid particlesP^ Similarly, the 
outer radius R may never reach the maximum allowed 
by the total arc length s if external conditions inhibit 
the scrolling. Below we consider two such conditions: 
mechanical clamping and electrostatic repulsion. 

When a CNS forms at the free edge of an initially 
flat strip, it is often clamped at two other edges, see 
Fig. [T] In this geometry, finite in-plane elastic deforma- 
tions build up as CNS rolls inward. Computing the cor- 
responding equilibrium shape of the CNS is much more 
complicated than for the uniform scroll, and so we make 
further approximations. We assume that everywhere ex- 
cept near the corners, s(x) and R(x) are slowly vary- 
ing functions of \x\ < L/2 that satisfy the inequalities 
s(x) ^> R(x) ^> r ~ sq. This allows us to treat the CNS 
as a slender curved beam. Its axis stays predominantly 
in the plane and its shape y(x) is dictated by the equilib- 
rium between the force —de/ds ~ 7 that pulls it in and 
the internal longitudinal tension T that tries to keep it 
straight: 7 ~ —Ty"(x). Thus, function y(x) is parabolic: 



y(x) ~ 2/(0) 1 
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For example, if L = 2 /im, then y(0) « 20 nm, in agree- 
ment with observations, see Fig. [T] 

The derivation of Eq. Q relies on the smallness of 
the ratio 7/Y2 <^ 1. While the unscrolled portion of 
graphene is stretched outward by the same force per 
unit length 7, the characteristic elastic deformation is 
u(x) ~ (j/Y 2 )L(l - 4x 2 / X 2 ) 1- " • Here the exponent 24 
a « 0.16 accounts for the power-law singularities near 
the corners p£l From Eq. (|3|, we see that u(x) is paramet- 
rically smaller than y{x), and so the unwrapped portion 
of graphene acts as an inextensible membrane. This im- 
plies that the "unstretched" (T = 0) shape of the beam 
is the parabolic arc located halfway between y = and 
the actual beam position. The tension T is related to the 
length difference AL of these two parabolas: 
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The denominator y(x) in the last integral takes into ac- 
count that the local cross-sectional area of the CNS is 
s(x)h with s(x) ^ y(x) for the inextensible graphene. 
The integral logarithmically diverges at the end points 
where our approximation is invalid. Imposing a short- 
distance cutoff d m i n , we recover Eq. (|3|. (The numerical 
answer given is nearly independent of this cutoff.) Fi- 
nally, the radius R(x) of the CNS can be obtained from 
s[x)\ R(x) ~ y/2as(x) ~ y/2ay(0)(l - 4x 2 /L 2 ). 

Let us now go back to the case of a uniform scroll but 
add interaction with substrate and external electrodes. 
Suppose the sample is deposited on an insulating sub- 
strate of thickness b and dielectric constant h. The other 
side of the substrate is covered with a metallic gate, see 
Fig.|3ja). If a voltage V is applied between the graphene 
and the gate, the total energy of the system acquires the 
electrostatic contribution e e \ = CV 2 /2 — qV = — CV 2 /2, 
where C and q = CV are the capacitance and charge per 
unit length of graphene, respectively. (The first term is 
the charging energy, the second one is the work of ex- 
ternal sources). Suppose that the CNS can form only at 
one edge and let I be the width of the unwrapped por- 
tion. Taking the fully wrapped state, I = 0, to be the 
energy reference point, we get 



s(l)=gl-[C(l)-C(0)]V 2 /2. 



(5) 



Since C(l) increases with Z, the charging of graphene 
would tend to unwrap the CNS. We will show that there 
is a critical voltage 



V, 



c2 
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(6) 



above which the unwrapping is complete, I = s. Here 
we assume that the CNS is thick enough, R ^> / 'g, 



3 



(a) 




K, V 



FIG. 3. (a) Model geometry (see text) (b) Optimal / vs. V. 



so that the effect of the bending modulus D on g can be 
neglected. However, g = 7 — 7 S can be reduced from 7 
due to the van der Waals coupling 7 S to the substrate. 

For a quick estimate, consider the standard Si02 sub- 
strate {h = 3.9, b = 300 nm). Unfortunately, reported 
7 S range widely, from about one half 26 27 (similar to 
metal substrates 28 ) to one 

hundredtrPSo] G f 7 itself. 
Presumably, 7 S is strongly influenced by surface condi- 
tions. Assuming 7 S <C 7, we get the upper bound on 
the critical voltage, V c2 = 65 V, which is still in the ex- 
perimental range. The corresponding carrier density is 
n c2 xV c2 /(^7reb) 4.6 x 10 12 cm 2 . 

In order to derive Eq. Q and find the optimal I at 
V < V c2 we need to have a model for the capacitance 
C(V). It is well known that C involves both classical and 
quantum terms. The former is the capacitance computed 
treating the CNS as a metallic cylinder whose charge 
resides only on its surface. Quantum effects make the 
charge to spread to a depth of the order of the screen- 
ing length, reducing the capacitance. This length can 
be estimated using the Thomas-Fermi (TF) approxima- 
tion. Locally, the CNS is similar to a stack of 2D lay- 
ers with the total density of states v\jh where v\ is the 
density of states of a single layer. The TF screening 
length of such a stack is ttf = yV (4ttz^i) = ^r\\hj2^ 
where r\\ = (2ne 2 vi)~ 1 is the in-plane screening length 
(for refinements, see Ref . I3T]) . For n ~ 10 12 cm -2 we have 
rn ~ 10 nm, so that ttf ~ 1 nm w 3h. Thus, for CNS 
of four or more coils, the classical capacitance is a good 
first approximation. Similarly, if I ^> ry, the unwrapped 
portion of graphene can be treated as a metallic strip. 

Assuming both such conditions are met, we can use 
some standard expressions for C(l). If I ^> 2b ^> i?, then 
C(l) is dominated by the flat piece and is given by the 
parallel-plate capacitor formula 
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(exact values of t^^tt ~ 1 are not important here), in 
which case Eq. ([5| takes the form 



4tt 2 



In 



const , 7 = 5 



xV 2 
8irb 



(8) 



In complete analogy to electrowettingP^the applied volt- 
age reduces the effective surface tension 7. The critical 



voltage V c2 [Eq. is determined simply by j(V c2 ) = 0. 
The origin of the force that unwraps the CNS is the asym- 
metry of the fringing field of the capacitor. This field is 
stronger on the outward side of the CNS because there 
it is not screened by the flat piece of graphene. Hence, a 
net torque on the CNS is produced. 

At V < V c2 the surface tension 7 is positive and the 
scrolling becomes possible. Minimizing e(Z), we obtain 
the optimal I as follows: 
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™ (v c2 /vy 



1 



(9) 



For k = 3.9, the range of validity of this formula is nar- 
row, 0.95 < V/V C 2 < 1, because at lower V we have 
I < 2b. Nevertheless, it does indicate that if we start 
with a flat sample of say micron dimensions at V = V c i 
and reduce V by just 5%, a rapid growth of the CNS at 
the expense of the unwrapped portion would occur. 

At still lower voltages, V < 0.8F C 2, at which I should be 
several times smaller than 26, we can use another analytic 
approximation: 



C(l) 



x+1 



4 In (8c(x)6/Z eff ) 
where l e g is the effective width of the strip plus CNS: 



(10) 
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ifZ»#, (11) 
if Z<ii. (12) 



We derived these equations using the standard tech- 
nique of conformal mapping followed by the series ex- 
pansion to the order indicated. In the limiting cases 
R = and I = 0, they reproduce the known results 
for the metallic strip 33 and the metallic cylinder, re- 
spectively. Function c{x) in Eq. (10) is given by a cer- 
tain cumbersome integral, which is however bounded and 
smoothly varying. Its characteristic values are c(l) = 1, 
c(3.9) = exp(-0.31) = 0.73, and c(oo) = 2/tt = 0.64. 
Analysis of the energy functional ([5| using Eqs. ( fToj ) 



( 12 ) unravels a new feature. There exists a lower critical 



voltage Vd, 
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/3 = 5.4, 



(13) 



at which the optimal width I discontinuously jumps from 
f3R to zero, i.e., the sheet becomes fully wrapped. For 
R 10-20 nm and b = 300 nm, we get V cl w 0AV c2 . 
Note that the I = state is a local energy minimum at 
any V. This is because at I <C R the field is symmetric 
on both sides of the CNS and no net torque exists to 
promote its unwrapping. This may also be the reason 
why no unwrapping has been seen in prior simulations 9 
even at carrier densities n ~ 10 14 cm 2 , i.e., hundred times 
higher than our n c2 . This leads us to believe that the 
transition to and from the fully wrapped CNS state would 



4 



be hysteretic in practice. Of course, no bistability would 
exist if the strip widths I < f3R are unattainable due to 
clamping or obstacles on the CNS path. Finally, Eq. (13) 



could be modified if the CNS length L varies as it wraps 
along, because of a non-rectangular sample shape. 

In closing, we posit that electrostatic and strain-based 
control over the shape and size of CNSs promise new ap- 
plications, such as "eyelid" actuators^ or valves in lab- 
on-a-chip devic es. Bein g ultralight, they would operate 
in a THz ranged 4 * 21 * 35 * Electronic paper and metamate- 
rials may be created with CNSs made of optically active 
multilayer graphene. Electronic applications may include 
var actors and THz oscillators p2 



Recently two other publications appeare d 35 * 36 * where 
the effect of longitudinal electric field on CNS was consid- 
ered. Similar to our case (the transverse field) it should 
create an accumulation of charge carriers, except they 



would have opposite polarity at the two ends of the CNS. 
The Coulomb repulsion of these mobile charges is an im- 
portant energy correction, which can be computed by a 
method similar to ours. In Ref . l36l only the polarization 
of immobile core states was included. We thank N. Pugno 
for drawing our attention to these papers. 
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